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In this work, we report electronic structure calculations aimed at computing the linear optical ab- 
sorption spectrum, and static dipole polariziblity of a newly proposed boron-based planar aromatic 
compound borozene {Bi2IIq). For the purpose, we use the semiempirical INDO model Hamilto- 
nian, accompanied by large-scale correlation calculations using the multi-reference singles-doubles 
configuration-interaction (MRSDCI) approach. We present detailed predictions about the energet- 
ics, polarization properties, and the nature of many-particle states contributing to various peaks 
in the linear absorption spectrum. Our results can be used to characterize this material in future 
optical absorption experiments. We also argue that one can deduce the aromaticity of the cluster 
from the optical absorption and static polarizability results. 

PACS numbers: 36.40.Vz, 31.10.+Z, 31.15.bu, 31.15.vq 



I. INTRODUCTION 



Next to carbon-based chemistry, boron chemistry is 
surely one of the most active areas of research^. Boron 
exhibits a variety of polymorphisms which involve three- 
center two-electron (3c-2e) bonds^, planar aromatic 
structures^, nano-tubes^, nano-ribbons^ and fullerene- 
like cages^. From the point of view of nano-technology 
aromatic boron materials are of particular interest be- 
cause delocalized electron clouds in these materials will 
lead to large linear and non-linear optical susceptibili- 
ties, just like corresponding carbon-based structures. In 
order to achieve aromaticity, one needs to have a planar 
structure giving rise to 7r-orbitals. Indeed, aromaticity 
in boron-based clusters has been the subject of many 
recent investigations. Zhai et al.'^ in a recent joint the- 
oretical and experiment study probed the issue of pla- 
narity and aromaticity boron clusters. In another theory- 
experiment study, Zhai et al.^ reported the discovery of 
planar boron clusters Bg and Bg clusters with "molecular 
wheel" like structures, and exhibiting both <j and tt aro- 
maticity. Aihara and coworkers^ presented a theoretical 
investigation of aromaticity in several planar and quasi- 
planar boron clusters using the concept of "topological 
resonance energy". Johansson^^ demonstrated theoret- 
ically the existence of strong magnetically-induced ring 
currents B20 and other toroidal clusters of boron, as an 
evidence of underlying aromaticity. Rincon et alr^ the- 
oretically demonstrated the a aromaticity in several pla- 
nar boron clusters performing an analysis of the electron- 
localization function (ELF). Similarly, Wu et al/^ using 
an ab initio approach, investigated the aromaticity of 
planar boron-carbon complex C6Bi2"^. Another possible 
approach to achieve planarity, and, therefore, aromaticity 
in boron clusters is by hydrogenation of boron clusters. 
For example Alexandrova et al}^ have demonstrated by 
means of theoretical calculations, that quasi-planar Bj 
cluster, upon hydrogenation, acquires a planar structure 
for B7H2 • In a recent theoretical work, Szwacki et alr^ 
proposed a novel planar aromatic cluster which can be 



obtained by hydrogenation of quasi-planar B12 cluster. 
They argued that the cluster is aromatic by performing a 
detailed analysis of its molecular orbit als, and by exam- 
ining magnetic properties such as nucleus-independent 
chemical shift (NICS), and anisotropy of magnetic sus- 
ceptibility (AMS). Having thus demonstrated the simi- 
larities between the aromaticities of this cluster and ben- 
zene (including the number of 7r-electrons) , they called it 
the boron analogue of benzene, and named it borozene^. 
Encouraged by the work of Szwacki et al}^ ^ Forte et 
al}^ performed ab initio calculations to predict larger 
aromatic compounds which can be constructed using 
borozene as formula unit. It is a well-known fact that 
the carbon-based aromatic materials including molecules 
such as benzene, napthalene, anthracene and longer tt- 
conjugated polymers exhibit intense linear and non-linear 
optical response^^. Therefore, it is of considerable inter- 
est to explore similar properties of aromatic materials 
consisting of elements other than carbon. With this aim 
in mind, here we present a systematic theoretical study, 
based upon a large-scale configuration interaction (CI) 
methodology, of linear optical response of borozene. We 
believe our results will be useful in optical characteriza- 
tion experiments on this material, as also in exploring 
its possible application in nano-optics. Furthermore, we 
also perform calculations of static polarizability of this 
cluster, with the aim of exploring the signatures of aro- 
maticity in its dielectric response. Aromaticiy is an intu- 
itive concept which essentially implies electron delocal- 
ization. Therefore, the consequences of this delocaliza- 
tion of electrons will be manifest in various properties of 
materials, including their magnetic as well as dielectric 
response. With the electrons delocalized along the plane 
of the molecule for conjugated systems as benzene, one 
would expect its "in-plane" dielectric response to be sig- 
nificantly more than its "perpendicular" response. Thus, 
this anisotropy in the static polarizability tensor aij can 
also be viewed as a signature of the aromaticity, in such 
systems. Indeed, this anisotropy in a has been verified 
theoretically by other authors in the past calculations on 
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conjugated molecules ^^^i^^i^^ and has been used to for- 
mulate an "aromaticity scale" for such systems.^'' 

Remainder of the paper is organized as follows. In 
section [J we briefly describe the theoretical methodology 
employed for the present calculations, and also discuss 
the geometry of the molecule. This is followed by the 
presentation and discussion of our results in section Hill 
Finally, in section [IV] we present our conclusions. 



II. THEORETICAL METHODOLOGY 

The basic structure of borozene is presented in Fig. 
[T] and it corresponds to point symmetry group Dsh—- 
Before proceeding with our calculations, for the purpose 
of veriflcation, we decided to perform the geometry opti- 
mization on our own using the B3LYP based hybrid DFT 
approach employing the 6-311++g(d) basis set as im- 
plemented in GAUSSIAN03 program^^. The optimized 
structure thus obtained belonged to the Dsh point group 
symmetry with uniform value of B-H bond length 1.18 
A, and four distinct values of B-B bond lengths 1.63 
A, 1.66 A, 1.81 A, and 1.86 A. We, however, are inter- 
ested in exploring an even more symmetric structure of 
borozene, so as to make the comparison with benzene 
more transparent. Therefore, we separately optimized 
a highly symmetric structure of borozene (S-borozene, 
henceforth) with identical B-B bond lengths, in addi- 
tion to the B-H ones. For this purpose, as well as for 
the linear absorption spectrum calculations, we adopted 
an INDO model Hamiltonian^^ based approach imple- 
mented recently in a computer program developed by 
INDO model is an effective valence-electron ap- 
proach, employing a Slater-type minimal basis set, and 
some semi-empirical parameters'^. The geometry opti- 
mization using our computer program, performed at the 
Hartree-Fock (HF) level (INDO-HF, henceforth), yielded 
bond lengths 1.65 A for the B-B bond, and 1.18 A for the 
B-H bond. These values agree perfectly with 1.18 A for 
the B-H bond length reported both by Szwacki et. ati^ 
and Forte et. al^^ , and the average value of B-B bond 
length 1.649 A reported by Forte et al^^ . The total 
energy of the geometry-optimized S-borozene was higher 
than that of the borozene with distinct B-B bond lengths 
by ^ 0.8 eV at the INDO-HF level. 

Next, large-scale correlation calculations were per- 
formed on the ground and the excited states of 
S-borozene, using the multi-reference singles-doubles 
conflguration-interaction (MRSDCI) approach as imple- 
mented in the MELD package'^, but employing the one- 
and two-electron matrix elements of the INDO Hamilto- 
nian supplied by our program^'. Thus, these correlated 
calculations were entirely within the INDO model, and 
win be called INDO-MRSDCI calculations, henceforth. 
During these calculations, point symmetry group C2v was 
utilized , as against the D^h symmetry group because the 
MELD'^ package is restricted to D2h and its sub-groups. 
Therefore, we classify the many-electron states of the sys- 




Figure 1: (Color online) Structure of S-borozene considered 
in this work, assumed to lying in the xy-plsbue. Optimized 
bond lengths at INDO-HF level are B-B = 1.65 A and B-H 
= 1.18 A. Yellow dots indicate boron atoms and violet ones 
hydrogen atoms. 



tem in terms of the irreducible representations (irreps) 
of C2v group. The optical absorption spectra was calcu- 
lated under the electric-dipole approximation employing 
the Lorenzian line shape. For both the ground and the 
excited states, these calculations were performed in an it- 
erative manner, until the optical absorption spectra com- 
puted from them converged. We have extensively used 
this approach in our earlier calculations on conjugated 
polymers^^ as well as on B12 icosahedral and quasi-planar 
clusters reported recently^^. 



HI. RESULTS AND DISCUSSION 

Symmetry group C2V consists of four irreps la- 
beled ai, a2, bi and 62- Assuming that the Carte- 
sian xy-pldine lies in the plane of the molecule (c/. 
Fig. [U, the symmetry adapted electronic structure 
of the HF ground state obtained from our code is 
la?2a?16?3a?26?4a?16p6?5a?46?6a?7a?8a?56?26|la|9a? 
6bl7b1l0al8bl, leading to ^Ai many-particle configura- 
tion. Some of the lowest unoccupied orbitals in the as- 
cending order of energy are 862, 2a2, 4^2, ll^i, 9bi etc. 
Of all the orbitals, the ones belonging to a2 and 62 irreps 
are tt orbitals and ai and bi are of a type. Thus we note 
that the highest occupied molecular orbital (HOMO) is a 
a orbital, while the lowest unoccupied molecular orbital 
(LUMO), is a TT orbital. For the sake of visualization, 
some of the orbitals are given in Fig. Illli 

While our LUMO orbital matches perfectly with that 
reported by Szwacki et. al^^, however, we have a dis- 
agreement on the nature of the HOMO orbital which was 
reported to be of 7r-type by them^^, but we obtain it to 
be of <j-type within the INDO model. In order to rule 
out the possibility that this disagreement could due to 
the use of the INDO model, or due to the S-borozene 
geometry, we studied the nature of the HOMO obtained 
using the 6-311++g(d) basis set both at RHF and DFT- 
B3LYP levels of ab initio theory^^ at various geometries. 




Figure 2: (Color online) Molecular orbitals (iso plots) of 
S-borozene from HOMO-1 to LUMO+1, obtained from the 
INDO-HF calculations. 



In all such DFT calculations, it turned out to be a cr- 
type orbital, but for the RHF calculations, the nature 
of HOMO was highly geometry sensitive with some ge- 
ometries yielding it to be <j-type, while others of 7r-type. 
This, in our opinion, is a significant difference between 
the electronic structures of benzene (whose HOMO and 
LUMO both are tt type), and borozene. Thus, borozene, 
in our opinion is an example of both a and tt conjuga- 
tion because, close to the Fermi level, occupied orbitals 
are mainly of <j type, while the unoccupied ones {LU MO, 
LUMO + 1, LUMO + 2) are of tt type. 

Before presenting the results of our calculations of lin- 
ear optical absorption spectrum of S-borozene, we would 
like to discuss the infiuence of geometry on the op- 
tical properties of borozene, considering the fact that 
our calculations were performed on a highly symmet- 
ric conformer S-borozene. For the purpose, we per- 
formed the optical absorption calculations at the singles- 
CI (SCI) calculations both on the lowest energy struc- 
ture of borozene with unequal B-B bond lengths, and S- 
borozene. We found that qualitatively, at the SCI level, 
spectra were very similar for both the geometries, except 
that for S-borozene it was slightly red-shifted compared 
to borozene. The first peak for borozene in this SCI cal- 
culation was obtained at 2.6 eV in perfect agreement with 
excitation energy of the first excited reported by Szwacki 
et alf^j while that for S-borozene was located at 2.4 eV. 

Next, we present and discuss the results of our INDO- 
MRSDCI based optical absorption calculations on S- 
borozene. Comparisons are also made with the absorp- 
tion spectra of quasi-planar B12 bare cluster, and also 
with that of benzene. Because the ground state of the 
system is ^Ai, as per electric-dipole selection rules, it can 
make linear optical transitions to excited states belonging 
to symmetry manifolds ^Ai, ^Bi, and through pho- 
tons polarized along the axis, x— axis, and axis, re- 
spectively. Thus, transitions to ^Ai and ^Bi are through 
photons polarized in the plane of the molecule (x?/-plane 
in the present case), while those to will be through 
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Figure 3: Linear optical absorption spectrum of S-borozene, 
computed using the INDO-MRSDCI approach. Important 
peaks, along with their polarization characteristics, are la- 
beled. A line width of 0.1 eV was used to compute the spec- 
trum. 

perpendicularly polarized photons. The INDO-MRSDCI 
calculations presented here involved construction and di- 
agonalization of very large CI matrices running into mil- 
lions of configurations, because simultaneously several 
states were targeted. For example, for the ^Ai symme- 
try manifold the total number of configurations was more 
than 2.1 millions, while those for ^Bi and ^B2 manifolds 
were in excess of 2.5, and 3 millions, respectively. Be- 
cause of the large-scale nature of these calculations, we 
believe that the results account for electron correlation 
effects properly. 

In Fig. [3] we present the combined optical absorp- 
tion spectrum of S-borozene, obtained from our INDO- 
MRSDCI calculations are presented. As is obvious from 
the figures that the first absorption feature is a small 
peak located at 2.25 eV, corresponding to a z-polarized 
transition into a B2 state. The many-particle wave func- 
tion of this peak consists mainly of two single excitations 
H ^ L^ldiudH-l^ L {H = HOMO; L = LUMO), 
both of which are of a ^ tt* type. After this feature at 
2.25 eV, the next peak follows after a large gap of ^ 5 eV, 
and is due to two closely placed states at 7.38 eV {Bi- 
type) and 7.55 eV (Ai-type). Thus, this intense peak can 
be reached through photons which are polarized in the 
molecular plane, but not through the z-polarized ones. 
The many-particle wave functions of both these states 
are dominated by singly-excited configurations, with the 
main configurations being \H^L + 3) for the Bi state, 
and \H — l^L + 3) for the Ai state. Thus, both these 
states correspond to <j ^ <j* transitions. 

Peak III of the spectrum, which is a weaker feature 
located at 8.44 eV, corresponds to a B2 state, and can 
be reached through a z-polarized photon. Many-particle 
wave function of this state is a mixture of several singly- 
excited configurations which include — 4 ^ L), \H — 
3 ^ L + 1), and ^ L + 8), all of which are a ^ tt* 
type excitations. Finally, feature IV near 9.1 eV draws 
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Figure 4: A comparison of the linear optical absorption spec- 
tra between S-borozene (solid line) quasi-planer B12 (broken 
line). The spectrum of S-borozene is significantly red-shifted 
as compared to that of B12. A line width of 0.1 eV was used 
to compute the spectrum. 



quasi-planer B12 cluster which has a convex shape with 
symmetry group Ca^^^, and can be seen as a precursor 
of borozene^^. Therefore, a comparison of the optical 
properties of quasi-planar B12 and borozene can help us 
understand the influence of hydrogen passivation on the 
optical properties of B12. In Fig. [H the linear optical 
absorption spectra of the two clusters is presented, and 
it is obvious that the flrst few peaks of S-borozene oc- 
cur at lower energies as compared to the quasi-planar 
B12. Moreover, the H ^ L transition is symmetry for- 
bidden in S-borozene, which is not the case for quasi- 
planar B12, and indeed contributes to the wave function 
of the flrst peak in its spectrum. Thus the main contribu- 
tion of hydrogen passivation on the optical properties is 
that it makes the structure completely planar (as against 
quasi-planar) leading to conjugation and higher symme- 
try. Lower energy gap of S-borozene is most certainly 
a consequence of conjugated nature of electrons in the 
system. 



its oscillator strength from the transition of the ground 
state to three closely spaced excited states of symmetries 
Ai, 5i, and B2 located at 9.03 eV, 9.01 eV, and 9.19 eV, 
respectively. The peak has mixed polarization features, 
out of which X- and polarized transitions are more 
intense as compared to the z-polarized one. The many- 
particle wave function of the Ai state is dominated by 
conflgurations \H — 1 ^ L + 5) and ^ L + 4), while 
that of Bi state consists mainly of — 1 ^ L + 4) and 
^ L + 5), ah of which are of a ^ a* type. The 
dominant conflgurations in the wave function of the B2 
state are \H ^ L -\- 8) and \H — 1 ^ L + 7), which are 
of (J ^ TT* type. The absorption spectrum of S-borozene 
consists of several very intense features beyond 10 eV as 
well, but we are not discussing them here, given the fact 
that the ionization potential of this material is estimated 
close to 9 eV— . 

Because Szwacki et alM called borozene the boron ana- 
logue of benzene, it is interesting to compare the optical 
absorption spectra of the two materials. It is a well- 
known fact that the flrst dipole allowed transition in ben- 
zene is from the ^Aig ground state to the ^Eiu excited 
state with the peak around 6.99 eV, through photons 
polarized in the plane of the molecule (x?/-plane in our 
case)^^. This compares very well with peak II (c/. Fig. 
[3]) of S-borozene which is also caused by photons polar- 
ized in the the xy-plane, and is located around 7.5 eV. 
The only difference being that for the the ^Eiu state in 
benzene corresponds to a tt ^ tt* transition, while the 
peak II of borozene corresponds to a ^ a* transitions. 
This supports the hypothesis that borozene is both tt and 
a conjugated^^. The difference between the optics of the 
two materials is the presence of low energy feature (peak 
I) in S-borozene corresponding to a z-polarized transi- 
tion, whose counterpart in benzene does not exist. 

Recently, we studied the linear optical absorption in 



In order to explore aromaticity of borozene, Szwacki 
et al}^ computed, and discussed its NICS plots, and 
anisotropy of magnetic susceptibility (AMS). With a sim- 
ilar aim we computed the three diagonal Cartesian com- 
ponents of the static polarizability tensor, axx^ c^yy, and 
azz of borozene both at the INDO-HF as well as INDO- 
SDCI levels. For the purpose we adopted a flnite-fleld 
approach, and computed the second derivative of the to- 
tal energy of the system numerically, with respect to the 
three Cartesian components of the external electric fleld, 
whose values was taken to be 0.005 atomic units. The 
numerical values obtained in our calculations correspond 
to valence electron contribution to the polarizability be- 
cause of the nature of the INDO Hamiltonian. However, 
our aim here is not to obtain exact numerical values of 
the polarizability tensor, rather to compare the values of 
its various components. It is obvious from the results of 
our calculations presented in table U that: (a) inclusion 
of correlation effects does not change the values of an 
very much, and (b) the values of in-plane components 
{axx and ayy) are signiflcantly larger than the perpen- 
dicular component azz- initio calculations performed 
on conjugated molecules like benzene ^^i^^ anthracenoi^, 
and ethane, exhibit similar anisotropy in the compo- 
nents of the polarizability tensor . And, indeed, based 
upon the magnitude of this "in-plane" vs. "perpendic- 
ular" anisotropy, Lazzeretti and Tossell^'' argued for an 
"aromaticity scale" to characterize 7r-conjugated systems. 
Using the same logic, the same anisotropy found in our 
calculations clearly implies the aromaticity in borozene, 
but as a consequence of both tt and a conjugation in the 
plane of the molecule. This fact is also obvious from the 
optical absorption spectrum of the material in which in- 
tensities of the z-polarized peaks are much smaller than 
those of X- and ^/-polarized ones. 
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Table I: Values of the components of static dipole polariz- 
ability tensor of borozene computed using the INDO-HF and 
INDO-SDCI methods (in atomic units). 



Method 


(^xx — ^yy 




INDO-HF 


34.273 


21.833 


INDO-SDCI 


33.705 


22.119 



IV. CONCLUSIONS 

In conclusion, a correlated study of electronic struc- 
ture and optical properties of a recently proposed boron- 
hydrogen cluster, borozene, was presented using the 
semi-empirical MRSDCI-INDO method. Our calcula- 
tions of the linear optical absorption spectrum can be 
used to characterize this substance in future optical ab- 
sorption experiments. We also calculated the diagonal 
components of the static polarizability tensor aij of the 



cluster, and found significant larger values for the com- 
ponents along the plane of the molecule, as compared to 
the one perpendicular to it. This anisotropy in the com- 
ponents of the polarizability tensor implies that the elec- 
trons are quite delocalized along the plane of the molecule 
leading to a considerably larger response to an electric 
field directed along the plane, as compared to that per- 
pendicular to the plane. Therefore, we argue that this 
anisotropy in the static dielectric polarizability of the 
cluster is an evidence of its aromaticity,^'' caused both 
by a and tt conjugation along the plane of the molecule. 
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